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We address the question as to why, in the semiclassical limit, classically chaotic 
systems generically exhibit universal quantum spectral statistics coincident with 
those of Random Matrix Theory. To do so, we use a semiclassical resummation 
formalism that explicitly preserves the unitarity of the quantum time evolution 
by incorporating duality relations between short and long classical orbits. This 
allows us to obtain both the non-oscillatory and the oscillatory contributions to 
spectral correlation functions within a unified framework, thus overcoming a signif- 
icant problem in previous approaches. In addition, our results extend beyond the 
universal regime to describe the system-specific approach to the semiclassical limit. 
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1. Introduction 

One of the central ideas in the field of Quantum Chaos is that the quantum energy 
spectra of classically chaotic systems exhibit universal fluctuation statistics (Bohi- 
gas et al. 1984, McDonald & Kaufman 1979, Casati et al. 1980, Berry 1987): once 
their mean energy level densities are scaled to be equal, the spectral statistics of 
different systems coincide, forming a small number of symmetry-dependent univer- 
sality classes. This idea is supported by overwhelming experimental and numerical 
evidence for systems as diverse as atomic nuclei, Rydberg atoms and quantum bil- 
liards (Stockmann 1999). Theoretically, it emerges in the semiclassical limit, H — ► 0, 
when spectral correlations are measured on the scale of the mean energy level sep- 
aration. 

The fundamental reasons for universality are still not fully understood. Nev- 
ertheless, a phenomenological description of the specific universal forms taken by 
spectral statistics has been developed based on the following insight: accepting 
universality, one may consider ensemble averages rather than individual chaotic 
systems. This idea lies at the heart of random-matrix theory (RMT), developed by 
Wigner and Dyson to describe spectral statistics in atomic nuclei (Wigner 1959, 
Haake 2001). In RMT, one averages over ensembles of all matrix representations 
of Hamiltonians belonging to the same symmetry class, with, for convenience, a 
Gaussian weight. The pertinent ensemble for systems without any symmetries, com- 
prising of all complex Hermitian matrices, is referred to as the Gaussian Unitary 
Ensemble (GUE); for dynamics whose sole symmetry is time-reversal invariance the 
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Gaussian Orthogonal Ensemble (GOE) of real-symmetric matrices is the appropri- 
ate one. Ensemble averages then yield predictions for measures of spectral statistics 
such as the two-point correlation function of the level density of energy levels E n , 

p(e) = £„<s(e-e„): 

fl<f)s ?("( E+ ^M £ -^))- L <"> 

Here (.} denotes an ensemble average, replacing an average over the centre energy 
E for individual systems, and p is the mean level density. In the large-matrix limit, 
the random-matrix predictions involve power series of non-oscillatory and oscilla- 
tory contributions with coefficients c ni d n (Heusler et al. 2007) depending on the 
symmetry class: 

oo / 1 \ n 

i?(e)=Re^( C „ + d rl c 2 ") -j . (1.2) 

n=2 

For the GUE, c 2 = -d 2 = -1/2 and c„ = d n = for all n > 2. 

Establishing precisely when and why random matrix theory describes the spec- 
tral statistics of individual classically chaotic systems is one of the central problems 
in the field. To illustrate its difficulty, we remark that several fully chaotic systems 
are known, such as the cat maps (Keating 1991) and geodesic motion on com- 
pact arithmetic surfaces of constant negative curvature (Bogomolny et al. 1992), 
for which the spectral statistics do not coincide with any of the random-matrix 
forms; and there are examples of perturbed cat maps where the spectral statistics 
do coincide with one of the random-matrix forms, but not the form one would pre- 
dict based on the symmetries of the classical dynamics (Keating & Mezzadri 2000). 
In all of these cases the origin of the unexpected behaviour lies in the existence 
of arithmetical quantum symmetries that have no classical counterpart. There are 
currently no rigorously established necessary or sufficient conditions for the quan- 
tum spectrum of any individual classically chaotic system to exhibit random matrix 
statistics. For this reason, the relationship is usually conjectured to hold for generic 
systems, where the term 'generic' is sufficient to cover experimental applications, 
but as yet has not been given a precise mathematical definition. 

The fact that we appear to be far from a rigorous understanding of universality 
in spectral statistics motivates the development of heuristic approaches to justify- 
ing the random-matrix conjecture in the expectation that this will lead to further 
insights into the reasons underlying its success, and its limitations. This is our main 
purpose here. 

A natural starting point for understanding quantum consequences of classical 
chaos is provided by semiclassical techniques. (Other approaches are discussed in 
Haake (2001), Muzykantskii & Khmel'nitskii (1995) and Andrccv et al. (1996).) By 
stationary-phase approximation of Feynman's path integral, Gutzwiller (1990) was 
able to represent the level density as a sum over classical periodic orbits. The cor- 
relation function R(e) is thus a sum over pairs of periodic orbits a, b (Berry 1985). 
Importantly, this classical double sum incorporates quantum interference effects 
through a phase factor e^ Sa_Sb ^ fi depending on the difference between the two 
actions S a and Sb- In the semiclassical limit, as ft — > 0, constructive interference 
requires small action differences S a — Sb, at most of the order of ft. Construc- 
tively interfering orbit pairs were identified by Hannay & Ozorio de Almeida (1984) 
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and Berry (1985), Sieber & Richter (2001), Sieber (2002) and Muller et al. (2004, 
2005), and shown to yield non-oscillatory contributions to R(e): Hannay & Ozorio 
de Almeida (1984) and Berry (1985) showed, using a sum rule that follows from 
crgodicity, that 'diagonal' pairs of identical or mutually time-reversed orbits give 
a contribution quadratic in 1/e and thus determine c 2 ; Sieber and Richter (2001) 
derived the cubic term C3 using pairs of orbits that differ from each other only by 
their connections in an 'encounter' of two close orbit stretches; and, building on 
that insight, it was shown by Muller et al. (2004, 2005) that higher-order contribu- 
tions to i?(e) are due to an infinity of further families of pairs of orbits differing in 
arbitrarily many encounters of arbitrarily many close stretches. Summation of all 
contributions yields the full non-oscillatory power series Re Y^=2 c ™ (~) ■ ^ n these 
calculations the essential simplifying assumption is that the orbit pairs identified 
are the only ones that ultimately contribute. We shall make the same assumption 
in our approach. 

The extension to the oscillatory terms proved to be much more difficult; these 
terms are related to more subtle correlations between periodic orbit actions (Arga- 
man et al. 1993). However, the leading oscillatory contribution was derived in a 
modified semiclassical setting by Bogomoly & Keating (1996), and an expansion 
to all orders realized by Heusler et al. (2007). The approach of Heusler et al. was 
based on the following detour inspired by RMT (or more precisely a field theoretic 
implementation of RMT, the so-called nonlinear sigma model): the correlation func- 
tion i?(e) can be represented through derivatives of a generating function involving 
a ratio of four spectral determinants A(E) = det(_E — H) with four different en- 
ergy arguments. (Here H denotes the Hamiltonian.) In fact, such a representation 
can be realized in two equivalent ways. It turned out that with a Gutzwiller-type 
approximation for the generating function, one either recovers the non-oscillatory 
or the oscillatory terms, depending on which of the two representations is chosen. 
This is somewhat paradoxical, because both representations should in principle 
give the same result. In each representation a part of the result is missed because 
in Gutzwiller's trace formula the energies have to be taken with sufficiently large 
imaginary parts to guarantee convergence, and these imaginary parts can, for ex- 
ample, cause oscillatory factors to become exponentially small. Nevertheless, it was 
observed that adding both results one recovers the full correlation function pre- 
dicted by RMT. This obviously leaves the question as to why the two types of 
contributions are additive, and why either one by itself does not give the complete 
answer. 

We here justify additivity within the framework of an improved semiclassical 
approximation developed by Berry & Keating (1990), Keating (1992) and Berry & 
Keating (1992). Their resummation of the Gutzwiller-type contributions preserves 
the unitarity of the time evolution (in particular, the fact that energy eigenvalues 
are real) and in doing so reveals a remarkable relationship between contributions 
from long and short orbits. It is this relationship that allows non-oscillatory and 
oscillatory contributions to be described naturally within a unified approach. 

We shall first review the resummation technique used and illustrate its appli- 
cation to spectral statistics for an averaged product of two spectral determinants. 
We will then turn to the generating function and the correlation function R(e), and 
finally justify the additivity assumed in Heusler et al. (2007). To keep the presen- 
tation simple, we will focus on systems without time-reversal invariance, however 
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the ideas generalize straightforwardly. Importantly, the treatment of the leading- 
order contributions will be carried beyond the universal regime, and so describes 
the approach to the semiclassical limit. 

2. Resummation 

To derive a semiclassical approximation for the spectral determinant one may start 
from Gutzwillcr's formula for tv(E — H) -1 . The spectral determinant is then ob- 
tained as an exponential involving the smoothed number N(E) of energy levels 
below E and a sum over classical periodic orbits a, 

A(E+) cx cxp(| £ rf £'tr-^-L_) 

cx cxp ( - «^V(£+) -J2 F ae lSJE+)/h )- (2.1) 

a 

Here, E + denotes an energy with a small positive imaginary part (needed to ensure 
convergence), and S a is the classical action of the orbit a. The factor F a incorporates 
an amplitude depending on the stability of a and the Maslov phase. It is convenient 
to expand the exponentiated sum over orbits into a sum over unordered collections 
of periodic orbits, or pseudo-orbits A (Berry & Keating 1990) 

A(E+) oc e -™N(E+)J2F A (-l) nA c iSA( > E+ V h , (2.2) 

A 

where ua is the number of orbits inside A and Sa is the sum of their actions. F A 
is the product of the factors F a of the individual orbits (it also includes corrections 
to the simple sign factor (— l)™- 4 that appear if A contains several identical copies 
of a given orbit). For negative imaginary parts, complex conjugation yields 

A(E-) oc j*"( E ~) FX(-l) nA e- iSA< < E ~V h . (2.3) 

A 

Importantly (2.2) and (2.3) do not manifestly preserve the unitarity of the quan- 
tum dynamics: in the limit of vanishing imaginary parts A(E + ) and A(E~) should 
become real (given the reality of the energy levels) and identical to each other. This 
is not at all obvious from the above formulas. It was shown in Berry & Keating 
(1990), Keating (1992) and Berry & Keating (1992) that by explicitly incorpo- 
rating unitarity one arrives at an improved approximation for A(E). Postulating 
that (2.2) and (2.3) become identical for E + ,E~ — ► E implies duality relations 
between pseudo-orbits whose duration (i.e., sum of periods of contributing orbits) 
is larger than half of the Heisenberg time T H = 2-Khp and those shorter than T/y/2; 
the overall contribution of long orbits is the complex conjugate of the contribution 
arising from short orbits. (In the case of billiards the relationship follows directly 
from the fact that the spectral determinant is an even function of the momentum 
and the spectral counting function an odd function, see Keating & Sieber (1994).) 
This leads to the improved semiclassical approximation for A(E) as a sum over 
pseudo-orbits with durations shorter than Tn/2, 

A(E) oc e - t7r7r{E) F A {-l) nA e iSA/h + c.c. (2.4) 

A (T A <T H /2) 
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Equation (2.4) incorporates explicitly the relations between long and short orbits. 
Taken together with the orbit correlations discussed above this leads to a complete 
semiclassical picture of spectral fluctuations. It is known as the Riemann-Siegel 
lookalike formula because of its similarity to a corresponding expression for the 
Riemann zeta function (Berry & Keating 1990). 

We emphasize that resummation solves an important problem inherent in (2.2) 
and (2.3): whereas (2.2) and (2.3) diverge unless the energy has a sufficiently large 
imaginary part, the sum in (2.4) is finite and so converges when E is real. Essentially, 
it is the need to include an imaginary part in the energy in (2.2) and (2.3) that 
is responsible for the analysis of Heusler et al. (2007) missing a part of the RMT 
expression for the spectral statistics, depending on the representation used. With 
(2.4) we are now able to access the complete expression. 

3. Product of spectral determinants 

To illustrate how resummation can be used in the context of spectral statistics, we 
first consider the product of spectral determinants 

11(a) = (A (E + a) A (E - a)} . (3.1) 

We shall see that for individual chaotic systems 11(a) conforms in the semiclassical 
limit to the random-matrix prediction (Ketteman et al. 1997) 

n(a) = 11(0) sin(27rpa)/27r7ja . (3.2) 

To show this, we insert (2.4) into (3.1). 11(a) then turns into the following double 
sum over pseudo-orbits A, B 

11(a) oc (c-^{E+ a) J2 F A (-l) nA e iSA ( E+a V H 

A (T A <T H /2) 

x ^N(E-a) F* B {-l) nB e- lSBiE - a)/h ^ +C.C. . (3.3) 

B (T B <T H /2) 

In (. . .} the non-conjugate terms in the approximation (2.4) for A(E + a) are paired 
with the complex conjugate ones for A(E — a). The term "+c.c." is due to the op- 
posite combination. All other pairings involve highly oscillatory factors with phases 
oc N(E) and thus vanish after averaging over E. 

Still, the remaining phase factor involving (Sa (E + a) — Sb (E — a)) jh oscil- 
lates rapidly in the limit h~ — > and thus averages to zero for most pseudo-orbits. 
Systematic contributions to (3.3) therefore arise only from pairs with action differ- 
ences at most of the order of h \. In the spirit of the diagonal approximation, the 
dominating contributions can be expected to originate from identical pseudo-orbits 
A = B. If we restrict ourselves to such pairs, and expand N (E ± a) ~ N(E) ±pa, 
S(E±a)n S(E) ± T(E)a, we obtain 

11(a) oce" 2 ^ \F A \ 2 e 2iaTA/h + c.c. (3.4) 

A (T A <T H /2) 

t We neglect contributions where the S- and jV-dependent phase factors both oscillate rapidly 
but these oscillations mutually compensate, since there is no systematic mechanism giving rise to 
the pairs of pseudo-orbits needed. 
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It is now tempting to drop the upper limit Tfj/2 since it tends to oo in the 
semiclassical limit and, given the factor l-F^I 2 , it is no longer needed for convergence 
(for a with an arbitrarily small imaginary part). Indeed, dropping the upper bound 
is justified immediately if semiclassically op — > oo. The result can then be written 
as 

11(a) oc c-^T^^ia/h) + c.c. (3.5) 
where the dynamical zeta function £(s) is defined by 

C( S )=^|^| 2 (-l)"' 4 e- sT - (3.6) 

A 

or, equivalently, 

r 1 ( S )=^|F A | 2 e- sT -. (3.7) 

A 

Importantly, in chaotic systems £(s) has a simple zero at s = (Haake 2001). This 
is equivalent to the periodic orbit sum rule that follows from classical ergodicity 
(Hannay & Ozorio de Almeida 1984). 

For energy differences of the order of the mean level spacing the same result 
is obtained after a short calculation: if we incorporate the condition Ta < Tfr/2 
through a step-function 8 (x) = J^-i^ ^e xz / h (with small positive c) H(a) can 
be written as 

11(a) oc e -*«P»J- I — e'^^lJU^e^^+cc. 

= e -2»7rpa 1 f C+l °° ^l e ^pz^ ( Z - 2l 
27ri Jc-ioo z \ h 

If we close the contour on the left, the residue at z = reproduces (3.5) and the 
residue due to the zero of ((s) at s = is proportional to (2ia)^ 1 + c.c. = and 
thus vanishes. 

We note that equation (3.5) includes non-universal corrections to the random- 
matrix prediction (3.2), encoded in the analytic structure of ((s). It reduces to 
(3.2) in the semiclassical limit because £(s) oc s as s — > 0; that is, as a consequence 
of classical ergodicity. We have thus seen how oscillatory contributions to spectral 
statistics arise naturally from a resummed semiclassical approximation, even beyond 
the universal regime. 



+ c.c. 



(3.8) 



4. Generating function 

Let us now consider the generating function for R(e), 

Here the two energy differences a, f3 in the denominator are taken with small posi- 
tive imaginary parts. Then 
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where (||) indicates the limit a,(3,j, 6 —> ^Pp- Equation (4.2) is easily checked: 

taking two derivatives of Z leads to ^ E+a-H ^ T e-0-h times the original ratio of 
four determinants; the two traces ultimately lead to level densities whereas the ratio 
converges to unity in the limit (||). 

To obtain a resummed semiclassical approximation for Z we use the Ricmann- 
Siegel lookalike (2.4) for the two determinants in the numerator. Due to the im- 
portant role played by the imaginary parts in the denominator, no resummation is 
possible there, and we rather stick to the unresummed expressions 

A(E + )~ 1 oc e"^( £+ ) F A e tS ^ E+ ^ h (4.3) 



and 

A(E-)- 1 oc e-"^^ FXc- lSA(E ' )/h . (4.4) 



Specifically, no resummation of A -1 as a finite sum can reproduce the poles at the 
energy levels (a finite sum can, of course, reproduce the corresponding zeros of A 
itself). The choice of when we use A(E + ) or A(E~) is fixed by the positions of 
these poles. Collecting the pseudo-orbit sums from all four spectral determinants 
we then obtain two sums over quadruplets of pseudo-orbits A, B, C, D, 



Z = I e ™N(E+ a ) ^ F A e iSA( - E+a 



)/h 



xe - t KN(E~p) p* e -iS B (E-P)/h 
B 

xe -MrlV(£;+7) p c ( y _lYc e iSc(E+ 1 )/h 
C (T c <T H /2) 

xe iTrN(E-8) ^ F^_l^D e -iS D (E-5)/h\ 
D(To<T H /2) I 
+{l -^S,S^ - 7 } . (4.5) 

These sums arise from combining the non-conjugate terms in (2.4) for A(E + 7) 
with the complex conjugate ones for A(E — S) and vice versa. 
Interference between pseudo-orbits leads to a phase factor 

e i(SA(£+a)-S B (£-/3)+Sc(£+7)-S , D(£-'5))/?i (4 g) 

The pseudo-orbits will interfere constructively if the cumulative action of A and 
C nearly coincides with the cumulative action of B and D. The simplest way to 
realize this is to have the orbits in B and D identical to those in A and C (neglecting 
repetitions). The sum over such "diagonal" quadruplets of pseudo-orbits can be split 
into sums over the intersections A n B, A n D, C CI B, C PI D. If we expand S and 
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N as above these sums can be written as 

ol 7rp(a+/3- 7 -5) \F AnB |V TAnB ( a +0)/ fi 



c 



^nB,AnD,cnB,cnD 
<.t c ,t d <t h ii) 



x |F An z5| 2 (-l)™- 4nD e iTAnD(a+<5)/;i 
x |F Cn s| 2 (-l)" cnB e lTcnB( '' + ' 3)/n 

x |fcnD|V Tc ^^+ 5 )/ n + {7^-(5, <5^- 7 }- (4.7) 

In the same way as for the product of spectral determinants, we now drop the upper 
limits at Ty/2f. This is clearly justified if semiclassically ap — > oo etc., and there 
are no singularities prohibiting analytic continuation to a, (3, 7, 6 oc 1/p. Evaluating 
the sum in that case gives 

«-i(a + [3)/H)«-i( 1 + 5)/H) +il °> d n{ ^> 

Like (3.5), equation (4.8) is valid even beyond the universal regime; that is, it 
describes the system-specific approach to the semiclassical limit via the analytic 
structure of the classical zeta function ((s). Taking the semiclassical limit we find 



7 — 7~ffu — T~r: + \7 ^ S, 5 ^ -7} = Z 1 + Z 2 



which is exactly the random-matrix prediction for the GUE (c 2 = —1/2,^2 = 
1/2). Here the summands Z\ and Z 2 are responsible for the non-oscillatory and 
oscillatory contributions to the correlation function, corresponding to two different 
saddle points in RMT. We stress that Z 2 becomes accessible to a semiclassical 
treatment only through resummation. Differentiating (4.8) gives a formula identical 
to that due to Bogomolny & Keating (1996) for the nonuniversal corrections to the 
random-matrix expression for R(e). 



5. Relation to Heusler et al. 

In Heusler et al. (2007) no resummation was made and thus only the first summand 
Z\ could be obtained. Hence if one represents R(e) as in (4.2) only non-oscillatory 
contributions are found. To recover oscillatory contributions an alternative represen- 
tation of R(e) through Z was used, with the limit (||) substituted by (x) a, (i — > 
7, 6 — > — 27775- If one replaces Z by Z\ (and drops — |) this time only oscillatory 
contributions are obtained. It was claimed that summation of both results 

2tt 2 p 2 dadf3 KID 2n 2 p 2 dadf3 k*) 2 [ ' ' 

f Again, convergence is guaranteed when a, /3, 7 and <5 are given arbitrarily small imaginary 
parts. 
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yields the full correlation function. After comparison to (4.2) we see that this claim 
is equivalent to 

^1 (5 2) 

which is easily verified from Ricmann-Siegel resummation: the right-hand side of 
(5.2) differs from g a g^ \ ^ because 7 and 5 are interchanged and flipped in sign; on 
the left the same effect is reached by replacing (||) with (x). 



6. Conclusions and outlook 

By resummation of the generating function we have clarified the semiclassical ori- 
gin of oscillatory contributions to spectral correlators. This opens the door for a 
systematic use of generating functions in semiclassical theory with possible appli- 
cations in quantum transport or to higher-order correlation functions. Moreover we 
have extended the semiclassical treatment beyond the universal regime. Within the 
diagonal approximation, non-universal effects come into play through the dynami- 
cal zeta function. It would be interesting to see if this carries over to off-diagonal 
terms. 

We are grateful to A. Altland, P. Braun, F. Haake, S. Heusler, M. Sieber and B. Simons 
for helpful discussions. 
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